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The 3a Faddeev equation using 2a RGM kernel involves redundant components whose 

' contribution to the total wave function completely cancels out. We propose a practical 

^ I method to solve this Faddeev equation, by eliminating the admixture of such redundant 

. components. A complete equivalence between the present Faddeev approach and a varia- 

' tional approach using the translationally invariant harmonic-oscillator basis is numerically 

' shown with respect to the 3a bound state corresponding to the ground state of ^^C. 
^ ■ 

o : 

CN ■ A first issue for applying realistic quark-model baryon-baryon interactions to 

. few-baryon systems, such as the hypertriton, is to find a basic three-cluster eqna- 

I tion which is formulated by using microscopic two-cluster quark-exchange kernel of 

■ the resonating-group method (RGM). This issue is non-trivial, not only because the 
^ , quark-exchange kernel is non-local and energy dependent, but also because RGM 
^ I equations sometimes involve redundant components due to the effect of the anti- 

• symmetrization; i.e., the Pauli-forbidden states. A desirable feature of such a three- 

• i-H , cluster equation is that it can be solved in either or both of the variational approach 

^ I and the Faddeev formalism, yielding completely the same result. In our previous 

■ paper EP, which is referred to I hereafter, we have proposed a simple three-cluster 
equation, which is similar to the orthogonality condition model (OCM),li* but em- 
ploys the two-cluster RGM kernel as the interaction potential. The three-cluster 
Pauli-allowed space is constructed by the orthogonality of the total wave functions 
to the pairwise Pauli forbidden states. Although this definition of the three-cluster 
Pauli-allowed space is not exactly equivalent to the standard definition given by the 
three-cluster normalization kernel, this assumption is essential to find a complete 
equivalence between the proposed three-cluster equation and the Faddeev equation 
which employs a singularity-free T-matrix derived from the RGM kernel (the RGM 
T-matrix) . a) 

More explicitly, the Faddeev equation for systems composed of three identical 
bosons is expressed as 

Xip = GofSif , (1) 
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where S implies the real symmetric i^S = S) matrix for rearrangement of the three- 
types of the Jacobi-coordinates: (^^ + = [(123) + (123)^] = Sipa- Note that 
1 + iS is semi-positive definite. In Eq. (1), (p = (pa is the a-component of the total 
wave function ^ = (pa + (fp + (f'y = {1 + S)ip, Gq = Gq{E) = 1/{E — Hq + iO) is 
the 3-body free Green function for the negative energy E, and T = T^^^ {E, e) is 
essentially the two-cluster T-matrix derived from the RGM kernel (see Eqs. (2.8), 
(2.9) and (3.20) of I). This T satisfies the basic relationship (Eq. (2.23) of I) 

{u\[l + Gof] = [l + fGo]\u) = . (2) 

Here we assume only one Pauli forbidden state \u), for simplicity. 
In order to find a trivial solution of Eq. (1), we first solve 

s\un = T\un , (3) 

where \uf) = \u)\f) is a product of two functions corresponding to the two momen- 
tum Jacobi-coordinate vectors, k and q, respectively. Each function of |n) and |/) 
is assumed to be normalized as {u\u) = {f\f) = 1- In the 3d' system discussed in I, 
there appears no r = — 1 eigen- value, while in the 3a system we find two solutions 
with r = — 1, as is shown below. We find these by solving the eigen- value problem 

{u\S\un=r\n (4) 

with smaller dimensionality than Eq. (3). In fact, if we have a r = —1 solution [f^) 
for Eq. (4), the relationship (1 + S)^ = [1 + (123) + (123)2]2 = 3(1 + S) implies that 
{u\{l + Sy\uf'') = 0. Ifwe further use S-t = S, this leads to {{1+S)uf''\{l+S)uf'^) = 
0, namely, (1 -|- S)\uf'^) = for r = — 1. Such a solution with r = — 1 corresponds 
to the [21] symmetric component of the total wave function with respect to the 
permutations of the three bosons. Let us assume that we only have one r = — 1 
solution for simplicity of discussion. The Faddeev component constructed from this 
solution 

ipl = Go\ur) with T = -l (5) 

gives vanishing contribution to the total wave function ^, since Go is [3] symmetric. 
On the other hand, it is very easy to show that ip^ is a \ = 1 solution of Eq. (1): 

(^5 = GofS^l for r = -1 . (6) 

This is because of the very special nature of T, satisfying the basic relationship 
Eq. (2). 

Since our true solution of Eq. (1) with A = — 1 includes such a component as ip^, 
a general solution of Eq. (1) with A = — 1 is expressed as 

p = {p + Cpl , (7) 

where 1^ is a special solution of Eq. (1) with A = — 1. An optimum choice of the 
coefficient C in Eq. (7) is to make the value {p\l + S\ip) / {ip\ip) max;imum. If we 
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expand if by the complete orthonormalized basis of Eq. (3) (including the r = — 1 
solutions), we find that this condition is expressed as the orthogonality 

{ur\(p) =0 for r = -1 . (8) 

Namely, we only need to solve Eq. (1) under the condition of Eq. (8). 

The following procedure to eliminate the redundant components with r = — 1 
from the Faddeev equation is very similar to the technique used for finding a unique 
solution of RGM equations. We modify Eq. (1) as 



GoTS - Go\un . . .^A uf 



<p. (9) 



If we multiply from the left, {uf^lGoTS = {uf^l yields 

x{ur\cp) = o. (10) 

Therefor, the solution (p of Eq. (9) satisfies Eq. (8) so long as A 7^ 0. The fact that 
99q is a A = solution of Eq. (9) implies that the trivial A = 1 solution of Eq. (1) 
is pushed down to the A = eigen-state in Eq. (9). The solution of Eq. (9) also 
satisfies Eq. (1) when A 7^ 0, since the second term of Eq. (9) does not contribute 
because of Eq. (8). If we multiply Eq. (9) only with {u\ from the left and leave the 
degree of freedom for q, we obtain 

\{u\^) = -{u\s^) - i^\Go\un ^^^^^^^^^^^^ {ur\^) . (11) 

Again, the second term of the right hand side does not contribute if A 7^ 0, and we 
find 

{u\X + S\ip) = for A 7^ . (12) 
In particular, the A = 1 solution of Eq. (9) leads to the condition 

{u\l + S\(p) = — *■ {u\(pa + (Pi3 + (Pf) = , (13) 

which implies that our total wave function does not contain the pairwise redundant 
component {u\^) = 0. In summary, solving Eq. (9) automatically guarantees the 
solution 

Xip = GofSip, {ur\ip)=0, {u\X + S\ip)=0 for A 7^ . (14) 



Table I. \uf} SU3 states for 3 a system. N = N1+N2 is the total h.o. quanta of the whole system. 
Only non-negative A and fi are allowed. 



iVi 


N2 


(Am) 





N 


(7V0) 


2 


N-2 


(NO), (iV-2,1), (TV -4, 2) 
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Table II. [3] and [21] symmetric basis for \uf), classified by the SU3 basis with (Xfi). 



N 


iVi 


N2 


(Am) 


[3] 


[21] 











(0 0) 


(00) 




2 





2 


(2 0) 








2 





(2 0) 


(20) 


(20) 


4 





4 


(4 0) 








2 


2 


(4 0), (21), (02) 


(40), (02) 


(40), (21) 


6 





6 


(6 0) 








2 


4 


(6 0), (41), (22) 


(60)2, (22) 


(41) 



Next, let us consider the analytic solution of Eq. (4) in a particular case of the 
3a system. It is convenient to use the translationally invariant harmonic-oscillator 
(h.o.) basis used in I. The Pauli forbidden states of the 2a system consist of the 
h.o. states with A'^i = (for the relative angular momentum ^ = 0) and Ni = 2 
(for 1 = and 2), where A^i is the total h.o. quanta (2n + 1) for the variable k. 
The odd A'^i state is omitted from the very beginning, by assuming only even partial 
waves for the 2a relative motion. We have to solve the problem which SU3 states 
are classified to the [21] symmetry among the S't/s-coupled 2-particle h.o. states 
[C/(jViO)(^)^(A'20)(Q)](AM)a- For A''! = and 2, these SU3 states are constructed as in 
Table I. Since we know that all the SU3 state of the 3a system is Pauli forbidden 
for the total h.o. quanta = A^i + N2 < 6, it is sufficient to consider only the states 
given in Table This result is obtained by enumerating all the [3] symmetric basis 
states using the Moshinsky's method™ tl For example, for the two independent 
(20) states constructed from Ni = 2, N2 = and A^i = 0, N2 = 2 h.o. bases, 
only one (2£j) state is [3] symmetric, while the other (20) state belongs to the [21] 
symmetryO Similarly, only one [3] symmetric (40) state is made for = 4, and 
another [21] symmetric Pauli forbidden state is constructed for the total angular- 
momentum L = states of the 3a system. These [21] symmetric SU^ states are 
explicitly given by 



(^[21] (20) 



^[21](40) 



1 

f2 



U(20)aiQ) - U(20)aik) 



U(m)a{ci) 



(20)fl 



U(20){k)U(^2Q){Q) 



(40)a 



(15) 



Note that the (21) and (41) SU3 states are not possible for L = 0. This analysis 
shows that the two r = — 1 solutions of Eq. (4) for the L = 3a system is nothing 



* Here, again, we can explicitly construct the basis states using the theory of Double Gel'fand 
polynomials, . 

** A simple discussion using the tensor components of the 3-dimensional Jacobi-coordinate vec- 
tors k and q also leads to this conclusion. Namely, the 21-dimensional representation of the N = 2 
2-particle h.o. states, kakp (6 dim.), Qaq/B (6 dim.), fec,q/3 (9 dim.), is decomposed into the following 
SU3 irreducible representations; [3](20): kak/s + qaQis (6 dim.), [21](20): kak/s — Qaq/s (6 dim.), 
[21](20): kaQp + Qckf} (6 dim.), [111](01): eaM^pq-i (3 dim.). The [21](20) state in Eq. (15) cor- 
responds to the second SUi state, which is [2] symmetric with respect to the exchange of the first 
two particles, but is not totally [3] symmetric. 
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but the [21] (20) and [21] (40) SU^ states in Eq. (15). It is also apparent why we have 
no r = — 1 solution for the 3d' system. In this system, the Pauli forbidden state 
is only (O.s) state for 2d', and there exists no [21] symmetric Pauli forbidden state 
constructed for = 0. 

It should be noted that the existence of the [21] symmetric trivial solutions in 
the original Faddeev equation (1) is essential to eliminate the three-cluster Pauli 
forbidden states, which can not be trivially eliminated only through the pairwise 
orthogonality conditions with respect to the variable k. The mechanism to guaran- 
tee such a favorable result is furnished by the cooperative role with the exchange 
symmetry of the boson system and the partial elimination of the functional space 
f{g) corresponding to the other variable q. 

It would be legitimate to ask why such [21] symmetric components admixed 
to the Faddeev component ipa do not play an important role in the usual Faddeev 
equations. Suppose we have a situation (pp + ipj = Sipa = —fa- Then the usual 
Faddeev equation, ipa = GoTa^fp + (p-y), becomes (pa = —GoTaPa- If one uses the 
relationship GgTa = GaV^, this equation is reduced to {E — -ffo)¥'a = 0, and pa 
turns out to be the plane wave. This implies that we have no square-integrable 
trivial solutions. Even if we use Xfa = GoTa{(pj3 + (p-y), we obtain 



which implies that there are no square-integrable solutions below the 2-body thresh- 
old. (Note that we are interested in the situation A > 0.) Since the [21] components 
admixed to ipa do not contribute to anyway, there is no need to worry about such 
admixture. 

For a practical application of the basic equation Eq. (9) to the 3a system, the 
self-consistency procedure for determining the energy dependence of the exchange 
term eK in the allowed space, discussed in I, is very important. We therefore need 
to evaluate the expectation value e of the 2a Hamiltonian through 



where Hq is the 3-body free kinetic-energy operator and the Faddeev component 
p is normalized as 3{(p\l + S\p) = 1. Starting from some specific values of e and 
E, we solve Eq. (9) and find a negative 3-body energy E such that the eigen-value 
X{E) becomes 1. The normalized Faddeev component (p yields a new vale of e 
through Eq. (17). Since it is usually not equal to the starting value, we repeat the 
process by using the new value. This process of double iteration converges very fast 
if the starting values of e and E are properly chosen. For numerical calculation, 
we discretize the continuous variables k and q, using the Gauss-Legendre ni- and 
n2-point quadrature formula, respectively, for each of the three intervals of - 1 - 
3-6 fm~^. The small contribution from the intermediate integral over k beyond 
ko = G fm~^ in the 2a T-matrix calculation is also taken into account by using 
the Gauss-Legendre n3-point quadrature formula through the mapping k = ko + 




(16) 



e = ^E + ^{p\Hoil + SM 



(17) 
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tan(7r/4)(l + x)n The momentum region q = 6 fmT - oo is also discretized by 
the 71-3 point formula just as in the k discretization case. The values of ni-n2-ns 
are given in Table III. The partial-wave decomposition of the 2a RGM kernel is 
carried out numerically using the Gauss-Legendre 20-point quadrature formula. The 
incorporation of the Coulomb force to the Faddeev formalism is only possible by using 
the shielded Coulomb potential n(r) = {l/r)9{Rc — r), where 6{x) is the Heaviside 
step function. Here we use a rather small value Rc = 6 fm, in order to make 
the phase-shift calculation of the 2a system for higher partial-waves numerically 
stable. The Coulomb exchange kernel for this interaction is explicitly calculated. 
Since our Faddeev component 93 oscillates due to the orthogonality to the redundant 
components, we need to include the 2a partial waves at least up to ^ = 4. The 
convergence of the result is confirmed by including the partial waves up to £ < 
8. The numerical inaccuracy of the angular-momentum projection is examined by 
extending the 20-point quadrature formula to the 30-point formula, and we found 
that the error is less than 1 keV if numbers of discretization points for k and q are 
large enough. Furthermore, the modified spline interpolation technique developed 
by Glockle et al.i^ is employed for constructing the rearrangement matrix S. For the 
diagonalization of the large non-symmetric matrix, the Arnordi-Lanczos algorithm 
recently developed in the ARPACK subroutine packageB^ is very useful. For the 
effective 2-nucleon force for the 2a RGM kernel, we use the Volkov No. 2 force with 
m = 0.59, following the 3a RGM calculation by Fukushima and Kamimurai^. The 
h.o. constant for the a cluster is assumed to be = 0.275 fm^^. 

Table III shows the solution of the Faddeev equation for the 3a system, obtained 
by solving Eq. (9). We find that the convergence is satisfactory, if the relative 2a 
partial waves up to imax = 8 are taken into account. We compare the best result of 
the Faddeev calculation with the variational calculation employing the [3] symmet- 
ric translationally-invariant h.o. basis. The model is explained in I. The numbers 



* These ns points for k are not included for solving the Faddeev equation (9), since it causes a 
numerical inaccuracy for the interpolation. 

Table III. 3a Faddeev calculation for L = 0, without the Coulomb force (with the shielded Coulomb 
force with Rc = 6 fm~^). The discretization points of k and q are specified by ni-n2-n3 (see the 
text), ^max implies the maximum partial waves included, and rimax = (^/2 + l)(3ni)(3n2 + ns) 
the dimension of the diagonalization for the Faddeev equation (9). The effective 2-nucleon force 
for the 2a RGM kernel is Volkov No. 2 with m = 0.59. 



■^max ni-n2-n3 


^max 


6(2a) 


£(3a 


) 


C(04) 


4 15-10-5 


4,725 


8.5826 (9.4542) 


-11.2023 (- 


-5.9313) 


0.8262 (0.7905) 


20-10-10 


7,200 


8.5825 (9.4541) 


-11.2032 (- 


-5.9318) 


0.8262 (0.7905) 


6 15-10-5 


6,300 


8.4490 (9.3511) 


-11.4151 (- 


-6.1188) 


0.8212 (0.7861) 


20-10-10 


9,600 


8.4488 (9.3509) 


-11.4158 (- 


-6.1193) 


0.8212 (0.7861) 


8 15-10-5 


7,875 


8.4467 (9.3492) 


-11.4179 (- 


-6.1212) 


0.8211 (0.7860) 


20-10-10 


12,000 


8.4465 (9.3490) 


-11.4187 (- 


-6.1217) 


0.8211 (0.7860) 
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without the parentheses in Table IV indicate the results without the Coulomb force, 
while those in the parentheses when the shielded Coulomb force with Rc = 6 fm~^ 
is included. (Note that the 3a energy with the present shielded Coulomb force is 
slightly (about 150 keV) more attractive in comparison with the that of the full 
Coulomb calculation in Table I of I.) Table IV also lists the expectation value of the 
2a subsystem, e, and the overlap with the simple shell-model wave function C(o4) as 
well. Here the SU^ (04) wave function is expressed in the 3a cluster model as 



y,[3]{04) 



t^{40)(fc)t^{40)(9)J ^^^^^ = Y^R2o{k,bi)R2Qiq,b2)Y^oo)oik,q) 



4 3 
-^^Ri2{k,bi)Ri2{q,b2)Y(^22)o{f^,Q) + -Ro4{k,bi)Ro4{q,b2)Y^4A)o{k,q) , 

(18) 

where bi = 1/47 ^2 = 3/I67 with ^ = = 21/, Rni{x, z/), the radial part of the 
h.o. wave function, and Y(^x£^i{k,q) the coupled angular-momentum functions. The 
difference between the Faddeev calculation and the variational calculation is very 
small. In particular the difference of the 3a energies is less than 1 keV both in the 
Coulomb on and off cases. 

In summary, we have found that the Faddeev equation using 2-cluster RGM 
kernel, derived in the previous paper, El^ may involve redundant components, if multi 
Pauli-forbidden states exist for the 2-cluster relative motion as in the 3a system. For 
systems of three identical bosons, these components do not contribute to the total 
wave function, since they belong to [21] symmetry with respect to the permutations 
of the three bosons. We have proposed a simple method to eliminate these redun- 
dant components through diagonalization procedure of the rearrangement matrix 
in the 3-cluster model space with the 2-cluster forbidden components. The equiv- 
alence between this modified Faddeev equation and the variational approach using 
the translationally invariant harmonic-oscillator basis is numerically shown for the 
ground state of the 3a system. 

This work was supported by a Grant-in-Aid for Scientific Research from the 
Ministry of Education, Science, Sports and Culture (No. 12640265). 
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Table IV. The L = lowest eigen-values for 2a and 3a systems, obtained by diagonalization using 
[3] symmetric translationally-invariant h.o. basis. N stands for the maximum total h.o. quanta 
included in the calculation, and rimax the number of the basis states. The numbers without the 
parentheses indicate the results without the Coulomb force, while those in the parentheses when 
the shielded Coulomb force with Rc = 6 fm~^ is included. 



N 




E{2a) 




e 


£;(3a) 


C(04) 


4 


- 


20.7705 (23.1969) 














6 


- 


9.7474 (1.8262) 














8 


1 


5.0088 (6.9205) 


22.7245 


(25.1180) 


5.4467 (12.6290) 




1 


10 


3 


2.5254 (4.3239) 


15.9632 


(18.0165) 


-2.9705 (3.6021) 


0.9567 


(0.9528) 


12 


7 


1.1182 (2.8321) 


12.9965 


(14.8448) 


-6.6379 ( 


-0.3657) 


0.9160 


(0.9077) 


14 


12 


0.2649 (1.9122) 


11.3357 


(13.0258) 


-8.5469 ( 


-2.4866) 


0.8860 


(0.8736) 


16 


19 


-0.2787 (1.3140) 


10.3336 


(11.8894) 


-9.6227 ( 


-3.7228) 


0.850 (0.8490) 


18 


28 


-0.6396 (0.9073) 


9.6987 


;il.l403) 


-10.2698 


(-4.4949) 


0.8511 


(0.8315) 


20 


39 


-0.6676 (0.6203) 


9.2852 


;i0.6295) 


-10.6721 


(-4.9962) 


0.8415 


(0.8190) 


22 


52 


-1.0629 (0.4111) 


9.0119 


;i0.2742) 


-10.9281 


(-5.3310) 


0.8349 


(0.8101) 


24 


68 


-1.1899 (0.2544) 


8.8295 


;i0.0235) 


-11.0934 


(-5.5590) 


0.8305 


(0.8037) 


26 


86 


-1.2839 (0.1343) 


8.7071 


(9.8446) 


-11.2014 


(-5.7168) 


0.8275 


(0.7990) 


28 


107 


-1.3547 (0.0403) 


8.6246 


(9.7157) 


-11.2726 


(-5.8275) 


0.8255 


(0.7957) 


30 


131 


-1.4090 (-0.0346) 


8.5686 


(9.6222) 


-11.3198 


(-5.9061) 


0.8240 


(0.7932) 


32 


158 


-1.4511 (-0.0952) 


8.5306 


(9.5538) 


-11.3514 


(-5.9624) 


0.8232 


(0.7914) 


34 


188 


-1.4842 (-0.1448) 


8.5047 


(9.5034) 


-11.3727 


(-6.0032) 


0.8225 


(0.7901) 


36 


222 


-1.5105 (-0.1860) 


8.4869 


(9.4660) 


-11.3870 


(-6.0329) 


0.8221 


(0.7891) 


38 


259 


-1.5317 (-0.2205) 


8.4746 


(9.4381) 


-11.3968 


(-6.0548) 


0.8218 


(0.7884) 


40 


300 


-1.5487 (-0.2496) 


8.4661 


(9.4172) 


-11.4035 


(-6.0711) 


0.8216 


(0.7878) 


42 


345 


-1.5626 (-0.2745) 


8.4603 


(9.4014) 


-11.4081 


(-6.0831) 


0.8214 


(0.7874) 


44 


394 


-1.5741 (-0.2959) 


8.4562 


(9.3895) 


-11.4113 


(-6.0922) 


0.8213 


(0.7871) 


46 


447 


-1.5835 (-0.3144) 


8.4533 


(9.3804) 


-11.4135 


(-6.0991) 


0.8213 


(0.7869) 


48 


505 


-1.5914 (-0.3305) 


8.4513 


(9.3734) 


-11.4150 


(-6.1043) 


0.8212 


(0.7867) 


50 


567 


-1.5979 (-0.3446) 


8.4499 


(9.3680) 


-11.4161 


(-6.1082) 


0.8212 


(0.7865) 


52 


634 


-1.6034 (-0.3570) 


8.4489 


(9.3639) 


-11.4168 


(-6.1113) 


0.8211 


(0.7864) 


54 


706 


-1.6081 (-0.3680) 


8.4482 


(9.3607) 


-11.4174 


(-6.1136) 


0.8211 


(0.7863) 


56 


783 


-1.6120 (-0.3778) 


8.4476 


(9.3582) 


-11.4178 


(-6.1154) 


0.8211 


(0.7863) 
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-1.6153 (-0.3865) 


8.4473 


(9.3563) 


-11.4180 


(-6.1168) 


0.8211 


(0.7862) 


60 


953 


-1.6182 (-0.3942) 


8.4470 


(9.3547) 


-11.4182 


(-6.1179) 


0.8211 


(0.7862) 


62 


1,042 


-1.6206 (-0.4013) 


8.4468 


(9.3535) 


-11.4184 


(-6.1188) 


0.8211 


(0.7862) 


64 


1,145 


-1.6227 (-0.4076) 


8.4467 


(9.3526) 


-11.4185 


(-6.1194) 


0.8211 


(0.7861) 


66 


1,250 


-1.6245 (-0.4133) 


8.4466 


(9.3519) 


-11.4186 


(-6.1200) 


0.8211 


(0.7861) 


68 


1,361 


-1.6261 (-0.4185) 


8.4465 


(9.3513) 


-11.4186 


(-6.1204) 


0.8211 


(0.7861) 


70 


1.478 


-1.6274 (-0.4231) 


8.4464 


(9.3508) 


-11.4186 


(-6.1207) 


0.8211 


(0.7861) 
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8.4465 


(9.3490) 


-11.4187 


(-6.1217) 


0.8211 


(0.7860) 



